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t*"*- ' We consider a source of gravitational waves of frequency uj, located near the 

center of a massive galaxy of mass M and radius R, with to S> RT 1 . In the case 
1 of a perfect fluid galaxy, and of odd-parity waves, there is no direct matter-wave 

o: 

interaction and the propagation of the waves is affected by the galaxy only through 
the curvature of the spacetime background through which the waves propagate. We 
find that, in addition to the expected redshift of the radiation emerging from the 
galaxy, there is a small amount of backscatter, of order M/oj 2 R?. We show that there 
is no suppression of radiative power by the factor 1 + uj 2 M 2 /4 as has been recently 

•i-H . 

^ ■ predicted by Kundu. The origin of Kundu's suppression lies in the interpretation of 

a term in the expansion of the exterior field of the galaxy in inverse powers of radius. 
It is shown why that term is not related to the source strength or to the strength of 
the emerging radiation. 



Typeset Using REVTEX 



1 



I. INTRODUCTION 



If a source of gravitational radiation is located in or near a massive body the curvature 
of spacetime caused by that body may influence the generation and propagation of the 
radiation produced by the source. If, for example, a quadrupole oscillator is located at the 
center of a galaxy of mass M and radius R, we might guess that the effect of the surrounding 
galaxy on the radiation produced inside it is of the same order, M/R, as the characteristic 
Newtonian potential, at least in the case, M/R <C 1, that the galaxy is nearly Newtonian. 
(Here and throughout we use units in which G and c are unity.) In this paper we analyze 
just what these effects really are for high frequency waves. 

To do this we consider a fairly specific astrophysical configuration. We suppose that 
there is a source of gravitational radiation emitting waves at frequency u, and confined to a 
central source region r < r$ of the galaxy. We require that the source be small compared to 
the radius of the galaxy (r$ <C R), We also require the source frequency to be high enough 
so that there is a region, outside the source (r > r s ) which is deep inside the galaxy (r <C R) 
and in the wave zone of the source (r 3> u~ v ). In this region it is meaningful to talk about 
the gravitational wave flux well inside the galaxy. For a typical galaxy R ps 10 23 cm and 
M 10 16 cm, and for a kilohertz gravitational wave uo ps 10 _7 cm _1 . For these values the 
region we require is the range of radii satisfying max(r,g, 10 7 cm) Cr< 10 20 cm. 

For such a configuration the standard analysis tells us that as the high frequency waves 
propagate outward there are two effects of spacetime curvature that affect their passage. 
First the frequency of the waves is redshifted so that the frequency observed far outside 
the galaxy is reduced from that at the source roughly by the factor (1 + $ ), where $o, 
the central potential, is of order M/R. The second effect is associated with the meaning 
of the radial coordinate. If "r" is the usual (i.e., Schwarzschild) radial coordinate, then for 
waves radially propagating outward, the rate dr/dt is slightly less than unity, and there is an 
attendant gradual phase shift of the waves, of order uM, as viewed in the r coordinate. (For 
a discussion of propagation of gravitational waves, and the distinction between generation 
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and propagation for an "isolated" source, see Thorne 0). 

We consider what other effects influence the propagation of gravitational waves, and 
find that there are interactions between the spacetime curvature and the waves which are 
interesting as points of principle, if not of astrophysical importance. There is, however, 
a significant additional motivation for such a calculation, and a major motivation for this 
paper. Kundu f|,||] has recently argued that gravitational wave energy propagating out 
of a gravitational potential well will be reduced in intensity by the factor (1 + M 2 u 2 /4) . 
Because ujM can be large (of order 10 9 for the typical numbers given above), such a reduction 
of kilohertz gravitational wave signals originating in other galaxies would make detection of 
signals impossible and would be of crucial importance in connection with the detection of 
gravitational waves by instruments now being developed. 

The remainder of the paper is organized as follows: We start in Sec. II by outlining the 
mathematical origin of Kundu's argument that gravitational radiation is suppressed. We 
then describe the argument against suppression given by Kozameh, Newman, and Rovelli 

and its relationship to the present work. In Sec. Ill we derive the necessary connections 
between the NP formalism, used by Kundu, and the formalism of metric perturbations. 
We show, in the Schwarzschild exterior, how the NP projection \l/o is related to the the 
Zerilli function [f| in the case of even-parity perturbations, and in the odd-parity case to the 
function solving the Regge- Wheeler equation |J. For outgoing solutions of both parities, 
suppression factors arise in the relationship between the terms describing the strength of 
radiation, and the terms describing the apparent quadrupole moment. The subsequent 
analysis then takes advantage of the considerable simplicity possible in the odd-parity case. 
A model problem is defined with a central source of odd-parity waves which propagate 
outward through a perfect-fluid galaxy. In Sec. IV a Green function solution to the odd- 
parity gravitational wave problem is constructed which shows clearly the relationships among 
the source strength of the waves, the intensity of the outgoing radiation, and the various 
terms that can be identified as the quadrupole moment. Section V takes up the problem of 
the extent to which the galaxy is transparent to (odd-parity) radiation. Numerical results are 
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then presented which show that even for strong gravitational fields, the effect of gravitational 
potential wells on the propagation of high frequency radiation is negligible (except, of course, 
for the well known redshift effect). A summary and discussion of conclusions is given in 
Sec. VI. 

II. THE KUNDU SUPPRESSION AND THE KNR MODEL 

Kundu's arguments are framed in the Newman- Penrose @ (hereafter NP) formalism and 
are based on the Weyl projection ty in that formalism. For an outgoing solution \l/ takes 
the form 

*o = ^oKMK 5 + 0(r- 6 ) (1) 

where u is retarded time. Due to its r~ 5 fall off at large r, the quantity \l/ is n °t usually 
viewed as a direct measure of radiation intensity for outgoing solution, but rather as encoding 
information about the multipole moments of the source in the near zone (i.e., at distances 
from the source small compared to a wavelength). The shear 

a = a (u,9,^r~ 2 + O(r- 4 ) (2) 

is well accepted as carrying the information about gravitational wave energy density, specif- 
ically in the Bondi news function || da^jdu. 

Kundu considers linear perturbations about a Schwarzschild background of mass M and 
shows that there is a simple relationship between the quantity ip® that carries information 
about multipole moments, and the quantity <Jq that carries information about radiation. To 
express this relationship it is convenient to define the "despun" equivalents \l/o and a, of 
the spin-weight +2 quantities and a, by 

*o = (l/2)99*o o = (l/2)99a , (3) 
where, on spin-weight +2 quantities, 
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For a multipole mode, of index I, in terms of despun quantities, we find 

d 3 $ l (£ + 2)\ dd d 2 d 

dv? 4 {£-2)\ du du 2 ' 1 ' 

in which the bar over do, in the first term on the right, indicates complex conjugation. 
A useful feature of despun quantities is that their real and imaginary parts correspond 
respectively to even- and odd-parity modes, so that, for £ = 2, we have from the real and 
imaginary parts of (Eh 
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When time dependence e is assumed, the result becomes 



<7Q 



(8) 



"6(1 ±zuM/2) 

with the + signs applying for even parity perturbations, and the - signs for odd. 

Kundu interprets this equation as telling us that the radiation amplitude, for a 
quadrupole source, is reduced due to the mass of the Schwarzschild background by the factor 
(1 ± iuM/2)~ 1 , so that the radiation power flux (proportional to \da /du\ 2 ) is reduced by 
his suppression factor (1 + o; 2 M 2 /4) _1 . 

Kundu's arguments depend crucially on his interpretation of as the quadrupole mo- 
ment of the source (aside from multiplicative factors). There are two types of justification 
given by Kundu for this identification. First, he argues that this identification is valid in 
linearized theory JTIJ , and is valid in the full nonlinear theory [DJ for stationary spacetimes. 
Furthermore, in the time dependent case ip® nas the required transformation behavior for 
the quadrupole moment. As a second and distinct justification, Kundu considers a gravita- 
tional wave source in a massive galaxy and argues that the source integral for the quadrupole 



moment will be affected by the galaxy only to order M/R, and therefore the quadrupole 
moment will be negligibly different from that if the source were in flat spacetime. 

A model problem has recently been published which suggests that Kundu's suppression 
factor is a mathematical artifact, and not of physical importance. Kozameh, Newman and 
Rovelli || (hereafter KNR) use a simple model of a perturbative scalar field $, and show 
that the same apparent suppression applies to quadrupole radiation in this model as in the 
case of gravitational waves. They consider the equations for i = 2 scalar perturbations in 
a Schwarzschild background of mass M and define ip with $ = if)(r, u)Y2 m (6, 0), where r is 
the Schwarzschild radial coordinate, u is retarded time, and Y-zm is a standard (spin-weight 
0) t = 2 spherical harmonic. An outgoing solution will then have the form 

in which the ipk are functions only of u. For a stationary solution if> and ipi would van- 
ish, and ip2 would be the quadrupole moment. To emphasize that the situation may be 
more ambiguous for dynamical solutions in curved spacetime, KNR refer to ip2 a s the "field 
quadrupole," and designate it by Qf . The equations for the scalar field in the Schwarzschild 
background then show that ipo, the radiative part of if) is related to Qf by 

^ = "T l+^M/6 - (10) 

For a given time changing quadrupole, therefore, the radiated power is reduced by the 
factor (1 + a> 2 M 2 /36) _1 from what it would be in flat spacetime. KNR assume that this 
suppression is analogous to the suppression of gravitational radiation found by Kundu, and 
that the scalar example provides a simple model for understanding the suppression. KNR 
then proceed to use the scalar model to investigate the question whether the field quadrupole 
Q* is really what is usually considered the "quadrupole moment" of a gravitational wave 
source, that is, whether Q* is the same as the "source quadrupole moment" Q. 

To address this question they consider a very simple model: a source inside a massive 
spherical shell of radius R. Since the spacetime inside the shell is flat, if) has a simple 
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closed-form outgoing solution for r < R. In the interior of the shell Qf can easily be shown 
to be identical to the source quadrupole Q. The interior and exterior solutions are then 
matched at r = R with the condition that the scalar field is continuous across the shell. A 
consequence of this is to require that Qf is discontinuous across the shell; it increases across 
the shell by the factor 1 + iuM/6. This enhancement factor cancels the suppression factor 
and one concludes that the relationship between the scalar radiation, and the scalar source 
quadrupole (aside from negligible factors of order M/R) is the same as in flat spacetime. 

The KNR model is very suggestive of the root of the problem, that the (exterior) "field 
quadrupole" differs from the source quadrupole. But one might ask what details of the 
model could be changed to make the argument more convincing. Two details would seem 
to deserve the most attention. First, the model involves scalar fields and it is difficult to 
be certain that the lesson of scalar fields applies to gravitational perturbations. A second 
detail of the KNR model is more important. In the KNR model the assumption that tp is 
continuous is tantamount to assuming that the shell is transparent to scalar radiation. But 
for a shell which cannot absorb or reflect radiation, in a time invariant background, we know 
a priori that the radiation outside must be related to the source in the same way as the 
radiation inside. The matching condition, then, eliminates at the outset any possibility of a 
Kundu effect. A more convincing calculation would model the interaction of the waves and 
galaxy to allow for interactions, in particular for backscatter. 

In the following sections we attempt to fill in some of these details. We study gravitational 
waves produced by a central source and propagating outward through a "galaxy." We include 
all effects of interaction and show that a nearly Newtonian galaxy is indeed transparent to 
the propagation of waves. We also show explicitly why the r~ 2 term in the outgoing solution 
does correspond to the quadrupole moment near the source at the center of the galaxy, but 
not outside the galaxy. 

III. RADIATION SUPPRESSION IN THE ZERILLI 
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AND REGGE- WHEELER EQUATIONS 



Although Kundu's analysis of gravitational radiation is carried out in the NP formalism, 
it turns out to be convenient, as well as instructive, to look at the problem in terms of 
metric perturbations. We start by showing the relation of the suppression factor in the two 
formalisms. Both inside and outside the galaxy we take the form of the background metric 
to be 



ds 2 = -e u dt 2 + e x dr 2 + r 2 {d9 2 + sin 



(11) 



with v and A functions of r only. We define the radial variable r* by 



dr/dr* = e 



(u-X)/2 = e a(r) 



(12) 



and the retarded time u by 



u = t - r*. 



(13) 



In the Schwarzschild geometry, even-parity perturbations for a particular multipole mo- 
ment (with £ > 2) are conveniently described by the Zerilli || function Z^ which satisfies 
a simple potential type equation. 



Z (+) = V + Z (+) j 

dr 2 dt 2 I 



where, for 



y( + ) = (1 - 2M M 

r 2 (l + 3M/2r) 2 



6M 9M 2 9M 3 



(14) 



(15) 



(We use here the notation of Eq. (62), Sec. 24, of Chandrasekhar |12]1 .) From this Zerilli 
equation one infers that for outgoing radiation has the form 



ZW = 4 +) ( M ) + 4 +) (w)r- x + z ( 2 + \u)r- 2 + 



(16) 



and that, for 



S 



dz{ +) /du = 3z { +) 
dz { 2 +) /du = z[ +) - 3M4 +) (17) 

dzi +) /du = -Mz{ +) + (21/4)M 2 4 +) - 

These equations show that Zq and z^ vanish for stationary £ = 2 perturbations. For 
stationary solutions the z^ term is the first nonvanishing term, and is considered to carry 
information about the source quadrupole moment. For nonstationary solutions it is the 
function z^ 1 that carries information about gravitational radiation, since the gravitational 
wave power is proportional to {dz^ / du\ 2 . From Eqs. |l7|we can then deduce a relationship, 
d 2 z 2 /du 2 = 3(zq —Mdzy /du), between the radiation quantity Zq\ and the "quadrupole 
moment" z 2 ■ This relationship implies a suppression factor |1 — iujM\~ 2 , which is different 
from the factor found by Kundu. The difference arises from the difference in the choice of the 
"quadrupole moment" one infers from and from , and demonstrates the importance 
of that choice in the inference of suppression of radiation. 

For even-parity perturbations, the relationship between Re(\I / o) an d is given in 
Sec. 31, Eq. (352) of Chandrasekhar [12j as 

Re(^o) sin 2 9 f 2M\ _ 3 + 3M/r + 9M 2 /2r 2 + 9M 3 /4r 3 ^ (+) 



£,-3/3 ^ r ) r 3 (l + 3M/2r) 2 



+ 



15 1- 3M/r - 3M 2 /2r 2 



rdu r 2 (l + 3M/2r) 

The following should be noted about our use of that result here: (i) The Zerilli functions 
defined by different authors differ by multiplicative factors, but overall multiplicative fac- 
tors will not affect the frequency dependent suppression factor, (ii) The Weyl projection 
\l/o is invariant with respect to infinitesimal tetrad rotations and infinitesimal coordinate 
changes. We therefore need not be concerned, for example, that Chandrasekhar employs 
a nonstandard coordinate gauge, (iii) Chandrasekhar assumes azimuthal symmetry for \l/o 
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and angular dependence C^ 3 / 2 (^) esc 2 6, which is proportional to the spin- weight 2 spherical 
harmonic 2^m, for m — 0. For this case the angular functions are pure real, and the real 
and imaginary parts of \I/o describe, respectively, even- and odd-parity perturbations. 

When the outgoing form for in Eq. [16] is substituted in the righthand side of Eq. [TR 



and Eqs. [H] are used, we find 

_«^ = £M +0(A (19) 

^1+2 

Here the quadrupole moment q(u) is given by 



q(u) = z ( 2 + \u) + (3M/2)z[ + \u). (20) 

Aside from numerical multiplicative factors, q(u) is equal to i/jq, and is what Kundu interprets 
as the source quadrupole moment. We note that Eqs. [T7| and |20] give us d 2 q/du 2 = 3[4 +) + 
(M/2)dzQ + ^ /du], and hence the even-parity Kundu suppression factor (1 + iuM/2)^ 1 . 

The function Z^> (in the notation of Chandrasekhar) , which describes odd-parity metric 
perturbations in the Schwarzschild geometry, satisfies the potential type "Regge- Wheeler" 
H equation 

(£- £V» = *->*->, (2D 



where 



V (~) = (1 2M/r) w + gM/r] (22) 

For a solution of the form 

= z { - ] + z ( - ] /r + z { 2 ] /r 2 + z { f ] /r 3 + • • • , (23) 
the Regge- Wheeler equation tells us, for £ = 2 multipoles, that 

dz^/du = 3zq~^ 

dz ( i~ ] /du = z[~ ] - (3M/2)4 _) (24) 
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dzi ] /du = 0. 



The relationship of an d ^ \ f° r odd-parity perturbations, is given in Sec. 31, Eq. 



(345) of Chandrasekhar [12 



as 



Im(dV /du) sin 2 9 ( 2M 



r -3/2 



1 

2^3 



£(£ + 1)- 



6M 



+ 



j_ / _ 3M | ^ g N 
r 2 1 r 9m , 



or 



(25) 



When the expansion in 



for an outgoing £ = 2 mode, is put into (ESI) we find 



Im(tftf oMO sin 2 5 / 2M 



a 



-3/2 
^+2 



,(-) 



+ C(r" 



(26) 



so that, aside from multiplicative constants, z% ^ is equal to dtp^/du and is the derivative of 
what Kundu identifies as the quadrupole moment. From Eqs. (Ell) we have that d 2 z^ /du 2 = 



3[zq~^ — (M/2)dzQ~' /du] and hence the odd-parity Kundu suppression factor (1 — iuM/2). 

The mathematics of the odd-parity modes can be much simpler than that for even-parity 
modes since the former modes do not couple to the perturbations of a perfect fluid. The 
issue of suppression is the same for both parities, so we choose to take advantage of the 
opportunity for simplicity and we consider below only odd-parity perturbations. 

In the standard formalism for metric perturbations, odd-parity motions are described as 
deviations of the metric in ([TT]) • We follow here the notation of Thorne and Campolattaro 
13| , in which the Regge- Wheeler || gauge is used and azimuthal symmetry is assumed. For 



odd-parity perturbations of multipole index £, the only nonvanishing metric perturbations 
are, in this notation, 



5gt<t> = ho(r,t) sin OdPi (cos 9)/ 39 



(27) 



^9r<f> = hx(r,t) sin 9dPe ( cos 9)/ 89 



(28) 



where Pi indicates the Legendre polynomial of index £. In terms of the notation of Thorne 
and Campolattaro, the Chandrasekhar function is 
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ZH = e a h 1 /r . 



(29) 



The Schwarzschild perturbation function is, of course, only defined in the exterior 
vacuum of the galaxy, where e a = e u = e~ x = 1 — 2M/r, but ( p9|) allows us to extend the 
definition of Z^~> to the interior. From the odd-parity field equations given by Thorne and 
Campolattaro, Z^~\ in the interior and exterior, is found to obey the following generalization 
of (Ell) 



e»£(l+l)-3r^e 2a 
dr 



Z<-> = s. 



(30) 



v dr 2 dt 2 

Here the source term S is defined, in terms of the perturbations SRe^ and SR^ of the Ricci 
tensor, by 



2re a 



d_ f e u 5R 6 
dr \ r 2 



e u sm 2 9 d ( 5R 



T(f> 



89 \ sin^ 6 



d_ 
86 



1 d 
sin 9 89 



P e (cos9) 



(31) 



We apply ([3(]) to the following model. At the center of a massive galaxy, of mass M and 
radius R there is a source of odd-parity gravitational waves at frequency uo. The source is 
confined to the region for r less than some source radius r^. The source, e.g., an oscillating 
neutron star, must of course not consist of a perfect fluid, since odd-parity perturbations do 
not couple to the motions of perfect fluids. For this reason we take the matter of the galaxy 
to be a perfect fluid, so that there can be a clean separation between the generation and 
the propagation of the gravitational waves, a separation that is not possible for even-parity 
waves. 

It is worth noting here that the resulting mathematical formulation differs very little 
from that for a scalar field of the type considered by KNR. Let scalar field $ have a source 
density £ so that 



(32) 



In the spacetime of (0), for a multipole of index £, this equation reads 



drl dt 2 J [r ' r 2 



(icy 

e v l(l + 1) + e 2a r^f- 
dr 



(r$) = re^S 



(33) 
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The difference between the form of this equation for r$ and ( p0|) for is only in the details 
of terms of order Mj R. We will show, in the next section, that these terms affect detailed 
numerical results but, for a nearly Newtonian galaxy, cannot cause significant suppression 
of radiation. Other features of the scalar and the odd-parity problems are parallel. In 
particular, for both cases we can consider a compact central source (no radiation originating 
from the bulk of the galaxy) and the matching conditions at the surface of the galaxy are 
that the fields and their radial derivatives are continuous. 

IV. ANALYSIS OF OUTGOING WAVES 

To investigate the nature of outgoing solutions we take the time dependence of the source, 
and of to be e lult , and we write 

ZH = tfj( r )e iuju S = e iu)t S(r) . 

The equation for odd-parity waves then takes the form 

i/;" + (a' - 2iuoe- a )i)' - ^[e x e(£ + 1) - 3ra> 

= e- 2a e iujr *S{r), (34) 

where ' denotes differentiation. 

It is straightforward, in principle, to construct a Green function solution to Eq. (|34"D 
from two homogeneous solutions. We define a "central" solution, ip c , as the homogeneous 
solution which is well behaved at r — > 0, with the limit 

Ur) ^ r e+1 . (35) 

The second solution is taken to be the "wave" solution if) w defined by the condition that it 
represents outgoing waves at large radii. The mathematical condition on this asymptotically 
outgoing solution is 
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lMr)^?l + 0(l/wr). (36) 

We define W = W(ip c: ip w ) = ^' w ^ c — i/j' c ip w to be the Wronskian of these two solutions, and 
we note that W must have the form 

W = e 2iujr *e- a(r) /K (37) 

in which K is a constant. 

If the source is confined to the region inside some radius rs, then for r > r s the Green 
function solution takes the form 



rrs 

Mr) = Kij w (r) / i) c {r)S{r)e- a e- lwr * dr 
Jo 



= K^{r)I u , (38) 

For r s <C R, and in the long wavelength (ur s <C 1) limit, the source integral I w has the 
approximate value 



/ r e+1 S(r)dr. 
Jo 



Note that the absence of a conical singularity requires A — > at r — > 0, but v{r — 0) = u 
will in general be of order M/R. It should also be noted that, aside from multiplicative 
numerical constants, the integral above is the usual integral for the £ th multipole moment of 
the source. 

The function ip w (r) corresponds to the solution that is asymptotically outgoing, but, 
due to backscatter, at small radius (r <C R), it does not in general have the appearance 
of a locally outgoing solution. We define a locally outgoing solution by the high frequency 
expansion 

/ x , ,o(r) b(r) .c(r) 

C^o u 

Here 
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oo = ooe- v °/ 2 



is the blueshifted frequency in the central region of the galaxy. This frequency governs the 
wavelength (A = 2irc/u) for r <C R) in the central region and is therefore the appropriate 
parameter to simplify (|3"5D. By solving the homogeneous wave equation ([34]) to various 
orders in ujq we find, for example, that a(r) and b(r) must satisfy 



1 



e a e- Uo/2 



.oc 



(40) 



b> = - e Q e- Uo/2 
2 



a + a a 



(41) 



The functions a, b,c,d, . . . are determined only after the metric functions z/, A are speci- 
fied, but we can state some general conclusions. For rCi? the metric coefficients can be 
expanded in powers of r and, for a geometry nonsingular at r = 0, we have v 1 = and A' = 
at r = 0. As a result, the solutions for a,b, . . . take the form 



+ 1)! M 



2r(£-l)\ R 2 



(42) 



1 (£ + 2)! M 
^{£-2)\ + R3 



v \ r 



(43) 



1 (£ + 3)! + M 



i? fr\ R 

c_i— In — + c_i h 

r Vit/ r 



(44) 



48r 3 (£-3)! i? 4 

and so forth. Here the coefficients a^, b^, . . . are numerical constants aside from corrections 
of order Mj R. More precisely they are functions of the parameters of the interior geometry 
which have finite limits as M/R — > 0. 



In the original homogeneous equation, Eq. ([H|) with the source set to zero, there is 
symmetry with respect to t —>■ — t and complex conjugation. From this symmetry we get a 
second, ingoing, solution 



(45) 
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in which the bar denotes complex conjugation. 

The asymptotically outgoing solution ip w (r) must be some combination of ip out and ip in , 
which we write as 



(46) 



where the constants T and 1Z can be considered transmission and reflection coefficients. 
The value of |T| 2 — \1Z\ 2 is computed by considering the Wronskian of = e lu}U ipw and 
its complex conjugate, and is found to be equal to unity aside from small corrections. (The 
value of \T\ 2 - \K\ 2 can be made precisely unity by a small correction in the normalization 
of ip w .) 



For r <C R the solution ip w (r) in ( f46|) can be expanded in powers of r. For £ = 2 this 



gives 



Til 



i 



0(Mr/R 3 



O ( — In 



i? 3 v# 



(47) 



+fte 2 ^ r * <i 1 



o; 



+ 0(Mr/R 3 



+ 



3 ^(Mir 
- + O ( — In 



For r < i? then, the term in ?/> w that goes as r~ 2 is 



^0 



(48) 



Thus, aside from corrections which are small for high frequency (o;i? ^> 1) sources, the r 2 
term in ^ is 



(T + ft) (-3/^)1^ 



(49) 
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If, as expected, backscatter is insignificant, then \T\ pa 1 and \1Z\ <C 1 so that the r 2 term 
is approximately 



-3/(4)KI u . 



(50) 



The r~ 2 term therefore gives a direct measure of the quadrupole source integral. From 
Eqs. ( p6|) and (|38]) it follows that ip — » ifi^ as r — ► oo, and that the r~ 2 term in the deep 
interior also gives a measure of the intensity of the outgoing radiation. 

For the exterior solution, very different conclusions follow. Here it is possible to expand 
ip w in inverse powers of r as 



+ — 

oo 1 



1 + - 



M 



3 , M A M 2 
- + A 1 — + A 2 — + 



-— + B X — + 

qy£i tyd 



1 

UO 6 



+ 



(51) 



in which the numerical constants, A\ = 3/2, A 2 = 0, B\ = 0, etc., are easily evaluated from 
the Schwarzschild metric functions. In the exterior solution then the r~ 2 term in ip is 



;-3/c; 2 + 3iM/2co)KI w . 



(52) 



and is larger than the interior r~ 2 term by the (possibly large) factor (1 — iuoM/2). But 
now the r~ 2 coefficient no longer gives the intensity of the outgoing radiation, or a measure 
of the source integral. If the coefficient is used to denote (in the notation of KNR) the 
"field quadrupole," then it must be understood that this field quadrupole is larger, by the 
factor (1 — iuM/2), than the quadrupole moment which measures the source integral, which 
governs the intensity of outgoing radiation at infinity, or which governs the locally outgoing 
radiation deep inside the galaxy. 

It is clear mathematically why the field quadrupole and the physical quadrupole are so 
different: the ia(r)/uo term in ([IT]) lacks a term that goes as r~ 2 . From d^0|) we see that the 
presence of such a term would require a galaxy spacetime that tends to a singularity as r — > 0. 
The absence of such a term is why, for a high frequency source, the r~ 2 term for r R 
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can be given the same meaning — that of the quadrupole moment — as in a flat spacetime 
background. In fl5T|) the ia(r)/uj term does have a r~ 2 term. This is possible because the 
expansion in fl5T|) cannot be extended to small r. But to interpret as the quadrupole moment 
the r~ 2 term, in some expansion for ip, is justifiable only if that expansion can be extended 
to small radii. The field quadrupole is, then, a formal construct and its use as a physical 
quadrupole moment is the reason that gravitational radiation appears to be suppressed. 

This insight gives the answer to an interesting question. Let us denote the coefficient of 
the r~ 2 term as Qf, both in the deep interior and in the exterior. When Qf is computed 
in the exterior we find a different value than in the interior, and than we would find in 
the absence of the galaxy. How can the "source integral" for Qf, when it is computed in 
the exterior, have large non-Newtonian contributions from the galaxy, especially in the case 
that the galaxy is nearly Newtonian? The answer is that the "source integral," both in the 
exterior and in the deep interior, is the same. It is I w of (0). But the way in which this 
source integral enters into the expression for the coefficient of the r~ 2 term, and hence into 
the value inferred for Qf, is different in the exterior and the interior. 

In the next section we consider just what the magnitude is of the influence of the galaxy 
spacetime on the outgoing radiation. 

V. REFLECTION AND TRANSMISSION OF OUTGOING WAVES 

We take up here the question of what the actual influence is of the curved spacetime of 
the galaxy on the propagation of gravitational waves (specifically, of odd-parity gravitational 
waves). One obvious influence, of course, is the redshift which is built into the expressions 
for the radiation. The net power in terms of coordinate time t must be independent of the 
distance from the source. The locally measured proper time differs from coordinate time by 
e v l 2 and hence the locally measured power (proportional to the square of the time derivative 
of ijj) will differ from that far outside the source by the redshift factor e u . 

The question of other influences on the radiation is much less obvious, and there are 
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several ways in which it can be asked. One approach is to look at the relation of the 
outgoing radiation and the source as embodied in ( |31| ) and fl38|). This approach is most 
transparent if the source is taken to be compact, i.e., rs <C as well as rs <C R. In this 
case, the radial derivative of v will be smaller (by rs/R) than the radial derivative of the 
Ricci components, so we can approximate 

dr V 0n o*J r d9\ sin 9 
1 d 



sin 9 89 



-PA cost 



(53) 



Here the terms SR§^ and SR~^ are the perturbations of the Ricci tensor projected on an 
orthonormal tetrad, and are the quantities that would be computed (e.g., for a neutron 
star) by a nearby observer. We therefore write 

S = e 3 "°/ 2 <S locaJ (54) 

to indicate the relation of the source term referred to the coordinates of ( fTT| ) and the source 
term measured by a local observer. 

For the compact source, with corrections of order r^M/i? 2 and l/r$w ignored, fl38f) can 
be written 

V(r) = i^°I loca ^ w , (55) 

in which 

f r S 

^locai = / r l+l S Xocal dr (56) 
Jo 

is the source term that would be computed by a local observer. In the case of flat spacetime 
K defined by ( 571 ) is easily shown to be — (i) e LV e [(2£ + 1)"] _1 so that, finally, the relation of 
source and field can be written as 

V(r) = -(i)V[(2* + l^rXcaiV'we^con- (57) 
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with 



corr — 



(2e+l)\\(-i) e u- t K. 



(58) 



In this equation the influence of the galaxy is contained in the factor e u °n corr . 

There is another, rather different, way in which the influence of the galaxy can be viewed. 
One can ask what the relationship is between the outgoing radiation far outside the galaxy, 
and the radiation in the deep interior of the galaxy. In fl46f) this relationship is contained 
in the constants T and TZ, in which TZ describes, approximately, the fraction of the radia- 
tion reflected back towards the source, due to the galaxy's spacetime curvature. Roughly 
speaking, the magnitude of \TZ\, or of \T\ — 1, is a measure of the extent to which the galaxy 
is not perfectly transparent to gravitational radiation. It is only an approximate measure 
because there is, at the outset, a limit to the precision to which an observer can measure 
radiation as if in flat spacetime. The metric for a "flat" coordinate system over a region of 
size L will deviate from the Minkowski metric by corrections of order (L/R c ) 2 , where R c is 
the spacetime radius of curvature. For the galaxy spacetime R c ~ (ii^/M) 1 / 2 , so that over 
one wavelength there will be metric corrections of order M/uj 2 R?. One manifestation of this 
is that we have, from the Wronskian of ip w and e 2luJT *ipw, and the expressions in P7|) and 



The 0(M / uj 2 R?) correction factor is simple to compute, once v and A are specified, from 
the forms for a,b, . . .. The correction factor M/uj 2 R? will, in any case, be negligible (of order 
10~ 39 for kilohertz waves and ordinary galaxies). 

There is a close relationship between the two viewpoints above for looking at the influence 



(0), that 



T| 2 - \TZ\ 2 = 1 + 0(M/cu 2 R 3 ). 



(59) 



of the galaxy. The Wronskian in ([H]) can be written 



w(if> c , = rf(i, vw) + nw^ c , ^ in ), 



(60) 



and the Wronskians on the righthand side can be evaluated to give, for £ = 2, 
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W^cVv) = (T + 1l)e 2iujr *e- a l5 uo 2 e uo/2 \l + 0(M/lu 2 R 3 )] . (61) 



As in fl59|) , the 0(M/u 2 R 3 ) correction term is easily evaluated from fl39| ) and the small-radius 
forms of a(r), b(r), . . ., once the the metric functions v and A are specified. 
When we combine ( 161] ) with ([J7|) and (|5^|), for £ = 2, we have 



e~ 3uo/2 {T + (l + 0{M/u 2 R 3 )) . (62) 



The effect of the galaxy is then contained in two types of terms. There are terms of order 
M/uj 2 R 3 (e.g., in (^) and (^)) that are "local" in the sense that they can be computed 
from the small-radius solutions for ip. The second influence of the galaxy is through the 
coefficient TZ, and is not local. If there is any way in which a nearly Newtonian galaxy can 
have a significant influence on the propagation of high frequency waves, it is through the 
possibility that \1Z\ is not small. 

That possibility can, in fact, be ruled out with a WKB argument, but such an argument 
cannot easily tell us how small \1Z\ really is. To find this out we have numerically integrated 
the equation for ip starting in the exterior, at large r, with the expansion in fl4"?|). The 



integration to small radius was done with the method of Gear [Q, suitable to the stiff 
differential equation for ip. At the surface r = R, the field equations require that ip an d ip' 
be continuous. The 1Z coefficient was extracted from the numerically computed solution ip w , 



by using the flat spacetime solutions = 1 — 3i/ru — 3/r 2 u 2 , and ipf^ 1 = e 2i< ^ r *?/>out) and 
the computed quantity 



ftindex = W(^t, ^)/W(^t, < n at ) . (63) 



This can be evaluated at small r with the expansion in (|39|) to give 



index 



+ 2ulR 3 



( —\ ~ tT 
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(64) 



where a\ is the coefficient defined in (f42|). Note that (|6J) does not assume M -C R] it can 
be used for galaxies with relativistically strong gravity. 
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For a constant density interior (Schwarzschild interior) the value of a% is easily shown 
to be 6, independent of M/R, and we apply ( j64| ) to the results for 7^ in dex found with ( |63| ) 
from the numerically computed values of ip w . In Fig. 1 we show the real part of 7Z m d C x as 
a function of r, for the parameters M/R = 0.1 and ujR = 2000, and we compare it to the 
predicted expression in pO) for the best-fit values TZ = (0.95 + 4.32)10" 8 and T = e a2 \ 
Numerical runs with different parameters show that \1Z\ is proportional to M/u 2 R 3 (aside 
from higher order corrections in ujR), so that the TZ and T terms on the righthand side of 
(|64l) are of the same order. 

Figure 2 gives \TZ\ as a function of M/R, for different values of ujR. The plots clearly 
indicate that \TZ\ ~ ^jj^ as long as M/R <C 1 and uR ^> 1. When M/R is no longer small, 
it remains true that \TZ\ oc uo~ 2 R~ 2 , but the dependence on M/R must be read from the 
figure. When uj is not large compared to R~ l , the assumptions used in deriving fl64|) fail as 
does much of the meaning of "reflection." 

VI. SUMMARY AND CONCLUSIONS 

We have studied a configuration in which waves propagate on a spherical static curved 
spacetime background. In the geometric optics limit, the limit of infinite frequency, the 
only physical influence of the background is the familiar redshift of the waves. We have 
found effects for high, but finite frequency u, for propagation in a background of mass M 
and radius R. Most notably, we have found that the radiation reaching arbitrary distances 
is different from that emitted, by a fractional correction of order M/u 2 R?. The computed 
reflection coefficients in Sec. V may be considered the first corrections to the geometric optics 
limit. 

Some details of Sees. IV and V are specific to odd-parity gravitational waves, but with 
very minor modifications apply also to massless, minimally coupled, scalar fields. The 
generalization to even-parity gravitational waves is not immediate. For the odd-parity, 
or scalar, case the waves propagate on the curved background of the galaxy fluid, but 
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there is no direct matter-fluid interaction. For even-parity waves propagating through a 
perfect fluid galaxy, or for waves of either parity propagating through a region with more 
complex material properties, the matter will, in general, oscillate in response to the passage 
of the wave, and will retard and absorb radiation much as a dielectric material interacts 
with an electromagnetic waves. (The computation for even-parity waves through a perfect 
fluid galaxy would be relatively straightforward to carry out with formalisms in which fluid 



perturbations do not explicitly appear |15| , but the problem is made difficult by its four 
degrees of freedom.) These interactions, however, can be estimated reliably [JTHJ] and except 
for contrived circumstances will be very small. 

What do these results imply for the possibility of suppression of gravitational radiation, 
as predicted by Kundu The use of a model calculation specific to odd-parity waves is 

irrelevant. The suppression is inferred by Kundu from the external Schwarzschild geometry, 
in which the mathematics of even and odd-parity waves is essentially the same. Our analysis 
is rather specific to a particular configuration: a compact central wave source embedded in a 
massive spherical background. One might ask whether the suppression might apply to very 
different configurations, such as a source near a massive black hole. It would be strange, of 
course, if the suppression — inferred only from the Schwarzschild background — applied for 
one wave source and not another within that background. Barring that possibility, the gen- 
eral lessons of our spherical configurations should apply insofar as well defined questions can 
be asked about suppression. In particular, for a compact source, in the geometrical optics 
limit (wavelength-C all other length scales), the effect of background curvature should be 
only the standard redshift and the bending of the null geodesies (absent in the spherically 
symmetric case). For finite frequency we would expect the first corrections from the geo- 
metric optics limit to be of order (ujR c )~ 2 , where R c is the characteristic spacetime radius 
of curvature. 

The conclusions based on the configuration considered in this paper should then give 
the generally correct picture of the relationship of radiation and quadrupole moment. In 
that picture we are able to distinguish a number of different "quadrupole moments": (i) a 
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quadrupole moment given by an integral over the source (ii) the interior "field quadrupole" 
inferred from the coefficient of the r -5 term in the NP Weyl projection near the source 
(iii) the "field quadrupole" from the r~ 5 term far from the source of background curvature. 
We have shown that the quadrupole moments of types (i) and (ii) are the same (aside from 
corrections of order M/uj 2 R 3 ) and that they govern the outgoing radiation produced by the 
source, both deep within the galaxy and outside the galaxy. The exterior "field quadrupole," 
however, differs significantly from the other quadrupole moments. Its interpretation as a 
quadrupole moment is based on an expansion in inverse powers of r, and an identification of 
the expansion coefficients with those of similar expressions for r R. But the expansion in 
the exterior cannot be extended inward, so that the expansion coefficients do not have their 
usual physical meaning. In particular the "field quadrupole" is not the quadrupole moment 
of the source, and does not govern the radiation produced by the source. 
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FIGURES 

FIG. 1. 7^-indcx is plotted for the numerical runs (circles) and the value given by the analytical 
estimate eq. (68) (curve). The parameters of the run are R = 100, M = 10, to = 20. The error 
bars in the numerical points are of the order of the size of the circles. 

FIG. 2. The absolute value of TZmdex plotted as a function of M/R for two values of u>R. The 
circles are for loR = 1000 and the stars for loR = 2000. The continuous curves represent the 
approximation for small M/R given by \R\ JJ^a • 
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